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Abstract. The mechanism underlying the long-range interstitial migration of nonequilibrium impurity 
interstitial species was used to simulate arsenic redistribution in ion implantation. An excellent agreement 
of the calculated arsenic concentration profiles with experimental data allows one to assume that the 
migration of nonequilibrium arsenic interstitial atoms makes a significant contribution to the formation 
of a low concentration region on thermal arsenic diffusion. The arsenic concentration profile calculated 
for a temperature of 1050 Celsius degrees within the framework of this assumption agrees well with the 
experimental one. A number of parameters describing arsenic diffusion at 1050 and 1108 Celsius degrees 
have been obtained. 


1 Introduction 

Recently, modeling arsenic thermal diffusion has been carried out in pQ on the basis of the diffusion 
equation 
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Here C and C B are the concentrations of substitutionally dissolved arsenic atoms and background 
impurity of the opposite type of conductivity, respectively; x is the concentration of electrons , normalized 
to the intrinsic carrier concentration nf, h ((7, (7 s ) is the factor describing the influence of the built-in 
electric field on the drift of charged pairs; D Sx , Df 1 , and Df 2 are the partial intrinsic diffusivities of 
dopant atoms due to the interaction with neutral, singly, and doubly charged vacancies, respectively; Df x 
is the partial intrinsic diffusivity of arsenic atoms due to the interaction with neutral self-interstitials; fa 
and fa are the empirical constants that describe the relative contribution of singly and doubly charged 
intrinsic point defects to the impurity diffusion. 

The empirical constants fa and fa in expression © can be found from the best fit to experimental 
data. For example, in Ref. [2] a fitting routine based on different experimental data was used and the 
temperature dependences of the arsenic partial diffusivities Df, D~, D 2 ~due to neutral, singly, and 
doubly charged intrinsic point defects were evaluated. The dependences obtained in [2] allow one to 
calculate the values of fa and fa for different temperatures of thermal treatment. It was shown in [I] 
that the concentration dependence of arsenic diffusivity obtained in [2] provides excellent agreement with 
experimental data for impurity concentration close to the value of rii. On the other hand, if the arsenic 
concentration C » rii and the concentration profile of electrically active arsenic gets a “box-like” form 
(see Fig. 3 in mi a difference from the experimental data is observed in the local region where a strong 
decrease in the impurity concentration begins. Therefore, for solving the problems arising in modeling 
silicon doping with arsenic it seems more convenient to use the concentration dependence proposed in |3], 
which is characterized by neglecting the interaction of impurity atoms with doubly charged point defects. 
Indeed, in this case, good agreement with experiment is observed in the whole region of high impurity 
concentration and there is only a small difference in the region of low arsenic concentration [T|. It was 
supposed in [T| that this difference appears due to the direct migration of arsenic interstitial atoms. 

The goal of this work is to investigate the possibility of arsenic interstitial diffusion and to obtain 
complete agreement of modeling results with experimental data. 

2 Interstitial diffusion of arsenic atoms in ion-implanted layers 

The assumption about the interstitial diffusion of arsenic atoms during ion implantation was first made 
in (4j with the purpose to explain the experimental data obtained. In |4] the p-type Czochralski-grown 
silicon substrates of (111) orientation with a resistivity of 5 - 10 flcm were implanted with As 75 ions at 
120 keV for doses of lxlO 12 - lxlO 16 ions/cm 2 . The average dose rate was approximately 1 pA/cm 2 . 
The implantation was carried out at a room temperature with the beam at a 7° angle to the silicon 
surface in an attempt to minimize channeling. After implantation, the wafers were cleaned and annealed 
at 600 900 °C in nitrogen environment. The conductivity profile of the implanted layer was obtained 

by using the incremental sheet resistance technique. Irvin’s well-known data were used to convert the 
conductivity profile into the electrically active arsenic concentration distribution. It was shown by a spe- 
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cial experimental procedure that the “tail” is formed during implantation and not as a result of enhanced 
diffusion from the high concentration region during the early stages of annealing. After annealing at 
600 °C, there was no evidence of an electrically active arsenic “tail”. When the annealing temperature 
was increased up to 725 °C, the tail became electrically active. The lack of electrical activity at 600 °C 
indicates that the arsenic is not in its usual substitutional environment in the “tail” region. It means 
that vacancies are not the traps for migrating interstitial species. 

In Figs. [T] and [2] the arsenic concentration profiles measured in J3] after annealing at a temperature 
of 725 °C for 60 minutes are presented. The implantation doses are 5xl0 13 and lxlO 14 ions/cm 2 , 
respectively. The chosen values of the doses ensure the maximal impurity concentration below a solubility 
limit of arsenic in silicon that is equal to 9.85x10 s pm -3 for a temperature of 725 °C and below an 
equilibrium electron concentration n e = 9.32xl0 7 /im -3 for this temperature 0- In Figs. [□ and [2] we 
also present the results of modeling the interstitial migration of arsenic atoms during ion implantation. 
A model of interstitial migration and an analytical solution of diffusion equation proposed in [B] were 
used for simulation of “tail” formation. Taking into account the results of U, it was also supposed that 
the interstitial arsenic species is trapped by uniformly distributed sinks but these sinks are not vacancies. 
To provide a best fit to experimental data, the Robin’s boundary condition is imposed on the surface 
of a semiconductor. Use of this boundary condition allows one to describe the evaporation of arsenic 
interstitial species through the surface. 



Figure 1: Calculated concentration profile of electrically active arsenic after annealing at a temperature 
of 725 °C for 60 min. Implantation dose equals 5xl0 13 ions/cm 2 . Experimental data are taken from [4]. 
The dotted curve is the Gaussian distribution after implantation. 

The following values of the model parameters were used to provide the best fit of the calculated arsenic 
concentration profile to the experimental one. 

For implantation with the dose Q — 5xl0 13 ions/cm -2 (Fig. [T|). 

The parameters prescribing the initial distribution of implanted arsenic are: R p = 0.077 pm, A R p 
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Figure 2: Calculated concentration profile of electrically active arsenic after annealing at a temperature 
of 725 °C for 60 min. Implantation dose equals lxlO 14 ions/cm 2 . Experimental data are taken from [4]. 
The dotted curve is the Gaussian distribution after implantation. 

= 0.028 /(in. The parameters specifying the process of interstitial diffusion are: the average migration 
length of arsenic interstitial species Iai = 0.052 /im, the fraction of the arsenic atoms participating in 
the interstitial migration p AI = 50.1 %. 

For implantation with the dose Q — lxlO 14 ions/cm -2 (Fig. [2]). 

The parameters prescribing the initial distribution of implanted arsenic are: R p = 0.079 pm, A R p 
= 0.034 pm The parameters specifying the process of interstitial diffusion are: the average migration 
length of arsenic interstitial species Iai = 0.044 pm, the fraction of the arsenic atoms participating in 
the interstitial migration p AI = 70.4 %. 

It can be seen from Figs. |T| and [2] that there is excellent agreement between the calculated arsenic 
concentration profiles and the electrically active arsenic profiles measured in jJJ. The agreement of 
simulation results with experimental data indicates in favour of the long-range interstitial migration of 
nonequilibrium arsenic species. It is worth noting that the parameters prescribing the initial distribution 
of implanted boron are approximately equal to the values tabulated in [T|: R p = 0.0712 pm, A R p = 0.0248 
pm, Sk = 0.31, R m = 0.0677 pm. Here Sk and R m are respectively the skewness and the position of a 
maximum of impurity distribution as implanted which is described by the Pearson type IV distribution 
0 . 

3 Thermal arsenic diffusion 

It was shown in [1| that the values of intrinsic diffusivity Di = 2.55xl0~ 6 pm 2 /s and the parameters /3\ 
= 611 and /?2 = 44.94 calculated from the expression given in j2j for a temperature of 1108 °C provides 
an agreement with the experimental data {8j for impurity concentration close to the value of rii. In Ref. 
[8], thermal arsenic diffusion was carried out from a constant source on the surface of a semiconductor. 
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The arsenic concentration profile for 120 min thermal treatment measured in [5] by neutron activation 
analysis is presented in Fig. [3] On the other hand, the results of our simulation of arsenic diffusion 
based on the model of [3], which neglects doubly charged point defects, give a similar agreement with 
the experimental data [B]. The following parameters that describe arsenic diffusion were used in this 
simulation: D t = 2.507xl0 -6 ^m 2 /s, 0\ = 100, and 02 =0 0- This agreement is not surprising because 
the maximal value of \ at the diffusion temperature is equal to 1.78 and 0iX»02X 2 ■ Only for X = 
13.6 we have 0ix ~ p2X 2 - H means that for the diffusion process under consideration a contribution of 
doubly charged point defects is negligible. Then, due to the great values of 0\ in the models [ 5 ] and [ 3 ], 
the arsenic diffusivity can be approximated by the identical expression 

D( X ) =D l D c ( X )~D lX . (8) 

In Fig. [3] a similar numerical solution of the diffusion equation (JT]) for the doping process investigated 
in [8] is also presented. Using a variation of the parameters £>* and 0\ (02 = 0), we have achieved an ideal 
fit to the experimental data of [5]. It is worth noting that the obtained values of Di = 2.7xl0~ 6 /xm 2 /s 
and 0i = 4, especially/?!, differ from the values reported in [2] and [3]. Therefore, it is interesting to 
investigate the arsenic diffusion with a higher doping level. The results of simulation of high concentration 
arsenic diffusion are presented in Fig. U The arsenic concentration profile obtained in [S] for a diffusion 
temperature of 1050 °C and duration of 60 minutes is used for comparison. 



Depth (gm ) 

Figure 3: Arsenic concentration profiles formed by thermal diffusion from a constant source on the silicon 
surface. The solid curve represents the arsenic concentration profile obtained by numerical solution of 
the diffusion equation that takes account of the concentration dependence of arsenic diffusivity and drift 
of the charged species in the built-in electric field. Diffusion temperature is 1108 °C for 120 min. Filled 
circles are the experimental data of [8|. 

It can be seen from Fig. |4] that the arsenic concentration near the surface is approximately equal 
to 3x10 s p,m~ 3 that results in x ~ 18.4. Thus, arsenic diffusion is characterized by a strong nonlinear 
concentration dependence of effective diffusivity. On the other hand, arsenic concentration is lower than 
the maximal equilibrium electron concentration n e = 3.566 xlO 8 /rm~ 3 for the annealing temperature 
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Figure 4: Arsenic concentration profiles formed by thermal diffusion from a constant source on the silicon 
surface. The dotted curve represents the arsenic concentration profiles obtained by numerical solution 
of the diffusion equation based on the pair diffusion mechanism. The solid curve takes addition account 
of the arsenic interstitial migration. Diffusion temperature is 1050 °C for 60 min. Filled circles are the 
experimental data of [8]. 

under consideration [5]- Therefore, it is possible to assume that arsenic clustering does not play an 
important role in the retardation of diffusion. The values of simulation parameters obtained from the 
best fitting to the experimental concentration profile are: Di = 5.8xl0 -7 /im 2 /s and /?i = 4 (fa = 
0). The impurity concentration at the surface has been chosen equal to 3.2 xlO 8 It is worth 

noting that the obtained value of Di is practically the same as Di = 5.8134xl0~ 7 /im 2 /s calculated for 
a temperature of 1050 °C from the expression proposed in [2] . On the other hand, fa = 4 and fa = 0 
used in simulation differ significantly from the values reported in [2]: fa = 238.9; and fa = 12.7. It can 
be seen from Fig. 4 that the calculated arsenic concentration profile described by the dotted curve agrees 
well with experimental data, except for the region of a low impurity concentration. It is similar to the 
arsenic profile calculated in [T within a framework of the model of m with Di = 5.537x10 ' pm 2 /s, fa 
— 100, and fa = 0. 

Taking into account the results obtained, it is reasonable to neglect the contribution of doubly charged 
point defects in simulation of high concentration arsenic diffusion in silicon, because in this case the 
calculated profile describes more precisely the high concentration region and the region where a strong 
decrease in the impurity concentration begins [T] . To achieve the agreement with the experiment in the 
low concentration region of arsenic distribution, we suppose that generation of nonequilibrium interstitial 
arsenic atoms can occur during high concentration arsenic diffusion. Participating in the long-range 
interstitial migration, these atoms form a low concentration region in the arsenic profile. For example, 
such generation can occur in the region of the abrupt fall in the arsenic concentration due to the stresses 
arising in the layer between the highly doped and intrinsic silicon. It is worth noting that in this layer 
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intense dissociation of the “impurity atoms — intrinsic point defect” pairs also takes place, as follows 
from the mass action law. After dissociation, the arsenic atoms that previously formed pairs become 
substitutionally dissolved again. However, a small fraction of arsenic atoms can occupy interstitial position 
and participate in the long-range migration. To take into account this additional impurity flux, the 
equation for interstitial diffusion [6] can be used. It is worth noting that thermal diffusion of arsenic occurs 
at temperatures significantly greater than the temperature of ion implantation in the experiments of [3] . 
Therefore, in contrast to jl] , it is supposed in this paper that arsenic interstitials become substitutionally 
dissolved again via recombination with vacancies (the Frank-Turnbull diffusion mechanism [9]). Then, 
the system of equations describing arsenic diffusion both by means of the formation, migration, and 
dissociation of equilibrium “impurity atoms - intrinsic point defect” pairs and due to the long-range 
migration of nonequilibrium arsenic interstitials has the form 


dC 

dt 
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Here C AI is the concentration of nonequilibrium interstitial arsenic atoms; d AI and t ai are the diffu- 
sivity and average lifetime of these nonequilibrium interstitial atoms, respectively; G AI is the generation 
rate of arsenic interstitials per unit volume of the semiconductor. We use the stationary diffusion equation 
for interstitial arsenic atoms in view of their large average migration length l A i = V d AI r AI and small 
lifetime ( t ai « r p ), where t p is the duration of thermal treatment. 

A numerical solution of the system of diffusion equations © and (flUl) for the process investigated in 
[8] is presented in Fig. [4] by a solid curve. It can be seen from Fig. [4] that account of the migration of 
nonequilibrium arsenic interstitials allows one to explain the formation of the “tail” in the low concentra¬ 
tion region and to achieve complete agreement of the arsenic concentration profile with the experimental 
one. The values of the simulation parameters that provide the best fit to the experimental concentration 
profile are: Di = 6.2xIO -7 /xm 2 /s and fa = 4 (fa = 0). The impurity concentration at the surface 
Cs has been chosen equal to 3.1 xIO 8 /im~ 3 and the average migration length of nonequilibrium arsenic 
interstitials Iai = 0.013 /jm. It can be seen from Fig. [4] that the calculated arsenic concentration profile 
described by the solid curve agrees well with experimental data over the whole diffusion zone. It is worth 
noting that full agreement with the measured concentration profile has been achieved using the assump¬ 
tion of the long-range migration of nonequilibrium arsenic interstitials. This assumption is based on the 
results of the previous simulation of arsenic diffusion during ion implantation. On the other hand, there 
can be another way to describe the low concentration “tail” region of arsenic distribution. For example, it 
was shown in |10] that in the low concentration region a condition of a local thermodynamic equilibrium 
between the substitutionally dissolved impurity atoms, intrinsic point defects, and “impurity atom 
intrinsic point defect” pairs can be broken up. In this case, the nonequilibrium pairs will be playing the 
role of nonequilibrium arsenic interstitials. Moreover, it was shown in mi that a similar “tail” region 
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can be formed in the case of diffusion of equilibrium “impurity atom intrinsic point defect” pairs if 
the nonuniform distribution of intrinsic point defects in the neutral charge state arises. Such nonuniform 
distribution can be formed due to the intense dissociation of “impurity atoms — - intrinsic point defect” 
pairs in the region of abrupt fall of arsenic concentration. 

4 Conclusions 

On the basis of the mechanism of the long-range interstitial migration of nonequilibrium impurity intersti¬ 
tial species, simulation of arsenic concentration profiles measured in [4j after implantation and annealing 
of silicon substrates at a temperature of 725 °C for 60 minutes has been carried out. The arsenic concen¬ 
tration profiles calculated for ion implantation with an energy of 120 keV and doses of 5 x 10 13 and 1 x 10 14 
cm -2 agree well with experimental ones. Thus, the simulation results obtained confirm the assumption 
of 0] that the “tail” in the low concentration regions of arsenic profiles is formed due to the interstitial 
migration of arsenic species during ion implantation. The average migration length of arsenic interstitial 
species obtained from the best fit to the experimental concentration profiles are equal to 0.052 and 0.044 
/im for implantation with doses of 5xl0 13 and lxlO 14 cm -2 , respectively. 

Based on the results obtained, simulation of thermal arsenic diffusion has been carried out. It is shown 
that assumption of the long-range migration of nonequilibrium arsenic interstitial atoms allows one to 
explain the formation of the “tail” in the low concentration region for the case of high concentration 
arsenic diffusion. Complete agreement of arsenic concentration profile with the experimental one [8] has 
been achieved over the whole diffusion zone. The average migration length of nonequilibrium arsenic 
interstitials obtained from the best fit to the experimental concentration profile is equal to 0.013 fj, m. It 
is worth noting that the empirical parameter /3i that describes the relative contribution of singly charged 
intrinsic point defects to the impurity diffusion is equal to the same value for diffusion at temperatures 
of 1108 and 1050 °C. The value of f3\ = 4 obtained from the best fit to the experimental arsenic profiles 
differs significantly from the values used in [2] and [5]. 
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